Slice Stretching Effects for Maximal Slicing of a Schwarzschild Black Hole 
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Slice stretching effects such as slice sucking and slice wrapping arise when foliating the extended 
Schwarzschild spacetime with maximal slices. For arbitrary spatial coordinates these effects can 
be quantified in the context of boundary conditions where the lapse arises as a linear combination 
of odd and even lapse. Favorable boundary conditions are then derived which make the overall 
slice stretching occur late in numerical simulations. Allowing the lapse to become negative, this 
requirement leads to lapse functions which approach at late times the odd lapse corresponding to 
the static Schwarzschild metric. Demanding in addition that a numerically favorable lapse remains 
non-negative, as result the average of odd and even lapse is obtained. At late times the lapse with 
zero gradient at the puncture arising for the puncture evolution is precisely of this form. Finally, 
analytic arguments are given on how slice stretching effects can be avoided. Here the excision 
technique and the working mechanism of the shift function are studied in detail. 
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I. INTRODUCTION 



When a singularity avoiding lapse together with a van- 
ishing shift is used to evolve a spacetime containing a 
physical singularity, the foliation is of pathological na- 
ture A first effect referred to as "slice sucking" 
consists of the "outward" -drifting of coordinate locations 
as the corresponding Eulerian observers are falling to- 
ward the singularity. This differential infall leads to large 
proper distances in between neighboring observers, cre- 
ating large gradients in the radial metric function, the 
so-called "slice wrapping" . 

The overall effect is referred to as "slice stretching" 
and, being a geometric property of the slicing, is present 
independently of the existence of a numerical grid Q . Its 
appearance has unpleasant consequences when numeri- 
cally evolving say the Schwarzschild black hole spacetime: 
Due to slice sucking the coordinate location of e.g. the 
event horizon is found to move outward and the outer 
region continuously decreases in coordinate size. This 
not only is wasting numerical resources but might cause 
problems as outer boundary conditions assuming nearly 
flat space and implemented at a fixed coordinate location 
become inappropriate and fail. Furthermore, not being 
able to resolve with a finite number of grid points the de- 
veloping steep gradients in components of the 3-metric, 
numerical inaccuracies caused by slice wrapping force a 
finite differencing code to crash. 

Following up on earlier work done together with 
B. Briigmann, [5|,IH0> m the present paper one particu- 
lar singularity avoiding slicing is looked at, namely maxi- 
mal slicing corresponding to the condition that the mean 
extrinsic curvature of the slices vanishes at all times @- 
This geometrically motivated choice of the lapse func- 
tion has been used frequently in numerical relativity, for 
simulations of a single Schwarzschild black hole see e.g. 

naiiiiimmiilH. 



For Schwarzschild, and by including electric charge also 
for Reissner-Nordstrom, the maximal slices can be con- 
structed analytically H E E El 13 . For those space- 
times it is hence possible to examine slice stretching ef- 
fects on an analytic level. The discussion throughout this 
paper is restricted to Schwarzschild, but since the same 
notation as in |EIE0 * s use d, the results carry over in a 
straightforward way to Reissner-Nordstrom. 

In the following the maximal slices of the Schwarzschild 
spacetime are re-derived in the radial gauge and the 
transformation to an Eulerian line element characterized 
by a vanishing shift is given. Performing a late time 
analysis as in based on it is shown here that 
slice stretching arises at the throat of the Einstein-Rosen 
bridge [H. 

Assuming symmetry with respect to the throat in order 
to fix its location in Eulerian coordinates by an isome- 
try condition, it is then possible to quantify slice sucking 
and wrapping at the event horizon acting as a "marker" . 
Examples of spatial coordinates to be discussed in the 
context of even boundary conditions are logarithmic grid 
coordinates (explaining numerical observations of e.g. 
[ToL Ir^ P and isotropic grid coordinates (extending the 
study of 0). 

For boundary conditions other than the even ones, 
however, this analysis is more involved as one impor- 
tant example, the so-called puncture evolution, shows 0. 
Here black hole puncture data is evolved using a lapse 
with zero gradient at the "puncture" , i.e. the compacti- 
ficd left-hand infinity. The corresponding lapse function 
is referred to as "zgp" or puncture lapse. 

Focusing on boundary conditions where the lapse arises 
as superposition of odd and even lapse, two integrals 
characterizing the overall slice stretching arc introduced. 
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It is then shown that for "favorable" boundary condi- 
tions the slice stretching effects can occur arbitrarily late 
in numerical simulations. Here the lapse at late times 
has to approach the odd lapse, the latter corresponding 
to the static Schwarzschild metric and, being antisym- 
metric with respect to the throat, having negative values 
in the left-hand part of the spacetime. The numerical 
implementation of a lapse function which is partially neg- 
ative, however, has been found to be unstable in at least 
two examples |2(j, |21j • Demanding hence in addition that 
a "numerically favorable" lapse should be non-negative, 
for latest possible occurrence of slice stretching a lapse 
being the average of odd and even lapse is obtained. One 
should note here that the puncture lapse is at late times 
characterized precisely by this condition. Furthermore, 
a one-parameter family of boundary conditions "ranging 
from odd to even" is studied numerically and the analyt- 
ically predicted overall slice stretching is observed. 



II. MAXIMAL SLICES OF THE 
SCHWARZSCHILD SPACETIME 

A. Radial gauge 

Following 0| and starting from the Schwarzschild 
metric in Schwarzschild coordinates {t, r, 9, </>}, 



1 



2M 



ds z = -f(r)dt z + —-dr z +r z dn z with f(r) = 1 , 

f(r) r 

(1) 

the maximal slices of this spacetime are most easily de- 
rived as level sets of the form 



a = t — t(r, r) = const. 



(2) 



Here the hypersurfaces are labeled by time at infinity r 
and one has to examine the behavior of the normal 



Finally, analytic arguments are given on how slice 
stretching effects can be avoided: 



n = NV (t - t(r, r)) =N [ dt - ^dr 

or 



(3) 



By excising parts of the hypersurfaces containing the 
throat, the latter being the origin of slice stretching, it 
turns out that analytically the evolution of the metric 
outside the excision boundary can be expected to freeze 
in the limit of late times. Such a behavior has been ob- 
served in numerical simulations, see e.g. pil 12^]. 

Not implementing a "throat excision" technique j2^| 
but making use of the shift function, it turns out that a 
coordinate singularity in the conformal factor is essential 
in order to "hide" the diverging term of the overall slice 
stretching there. For the puncture evolution this working 
mechanism is studied in detail for a "model shift" freezing 
the evolution of the rescaled radial metric component and 
locking the right-hand event horizon. For logarithmic 
grid coordinates, however, the shift function can not be 
expected to cure slice stretching. 

The paper is organized as follows: In Sec.[H]the max- 
imal slices of the Schwarzschild spacetime are re-derived 
and the origin of slice stretching effects is pointed out. 
In Sec. I1III slice stretching effects are studied, concen- 
trating in Subsec. llII Al on even boundary conditions and 
discussing logarithmic and isotropic grid coordinates ex- 
plicitly. The overall slice stretching is quantified in Sub- 
sec. I1II Bl and boundary conditions are derived which 
make those effects occur late in numerical simulations. 
As examples the puncture lapse and a one-parameter 
family of boundary conditions are discussed. In Sec. IIVI 
analytic arguments are given on how slice stretching ef- 
fects can be avoided, analyzing in Subsec. TlV Al the tech- 
nique of throat excision, and discussing for shift functions 
in Subsec. IIVBI their failure and in Subsec. IIV CI their 
working mechanism for logarithmic and isothermal grid 
coordinates, respectively. The results are summarized in 
Sec. E| 



Making use of the line element the normalization N 
is fixed by demanding 



"'■""= w2 im +fir) (K 



= -1. 



(4) 



As pointed out in |18|, N can also be considered the 
boost function of the static Killing vector relative to 
a = const, 



N 



(5) 



Calculating the trace of the extrinsic curvature, K 
turns out to be given by 



K 



d 


-r 2 f(r) 


dt 


1 








dr 





(6) 



Demanding for maximal slicing K = 0, the term in the 
brackets of 10 obviously has to be a function of time only 
to be denoted by C(r). Hence 



-fr r) = C{T) 

dr f(rWr 4 f(r) + C 2 (r) 



(7) 



is found, which has to be integrated by imposing bound- 
ary conditions in order to obtain the level sets J3Jl. 

Furthermore, using the future normal of the foliation, 
Ua = -qV^t. and writing the static Killing vector as 



tained as 



=JVn" + f with 



ol(t, r) = N(t, r) 



M = 0, the lapse can be ob 
dt 



3t 



(8) 
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Here N by the normalization Q is given by 



N(r,r)=±\ f(r) + 



±- 



(9) 



when introducing for convenience the polynomial 

p c (r) = r 4 /(r) + C* 2 = r 4 - 2Mr 3 + C 2 , (10) 

the subscript C denoting its dependence on C and hence 
on r. Together with the radial metric component 



7(r,r) 



Pc(r-) 



(11) 



and the shift 



p{Tir) = -(r,rMr,r) c{T) {u) 



the maximal slices in the radial gauge 

ds 2 = (-a 2 + dr 2 + 2f3drdr + jdr 2 + rW (13) 

have been derived. 

Furthermore, as pointed out in |24|, the radial and 
angular components of the extrinsic curvature turn out 
to be given by 

K = -2% and Kt=K+ = %, (14) 



respectively. 



B. Eulerian observers 

Next the spatial Schwarzschild coordinate r on the 
maximal slices shall be substituted by a spatial coor- 
dinate z corresponding to Eulerian observers. Apply- 
ing a transformation of the form r = r(r, z), the lapse 
a = ol{t, r(r, z)) is still given by (JSJ) and the line element 
is characterized by a vanishing shift, 



ds 2 = -a 2 dr 2 + Gdz 2 



(15) 



In the context of maximal slicing, K = 0, by contracting 
the evolution equation for the extrinsic curvature, one 
immediately obtains the statement that for zero shift the 
determinant of the 3-metric has to be time-independent. 
Hence the singularity avoiding property of maximal slic- 
ing comes to ligh t as the variation of the local volume 
remains fixed |25j . For this reason one can make for the 
radial metric component the ansatz 



G(t,z) 



H{z) 
r 4 (r, z) 



(16) 



where the function H(z) depending on z only is deter- 
mined by the initial data. 



The coordinate transformation relating r and z is 
found by comparison of the radial part of \i'6\ and l|15|) . 
For fixed slice label C one can infer the ordinary differ- 
ential equation 



dr 



(17) 



C— const 



which can be integrated using the throat as lower inte- 
gration limit by 



V 4 dy 

Vpc(y) 



= ± / Vimdy. 



(18) 



Here the 



or "-" sign applies for the right- or left- 



hand side of the throat, respectively. Furthermore, rc 
and zc denote the location of the throat in terms of 
Schwarzschild and Eulerian spatial coordinates. 

Note that rc is found as root of the polynomial 
Pc{t), which implies C(r = 0) = when starting with 
the throat of the Einstein-Rosen bridge coinciding ini- 
tially with the event horizon at teh = 2M. The throat 
rc, describing the "innermost" two-sphere on a slice la- 
beled by C, never reaches the singularity at r = as 
r c, im = 3M/2 is found in the limit of late times with 
C approaching 



-V3M 2 



(19) 



as pointed out in [9| . Hence for the Schwarzschild space- 
time the singularity avoidance of maximal slices becomes 
apparent, c.f. corollary 3.3 of pfij . 

The coordinate location of the throat in Eulerian coor- 
dinates depends in general also on C and is hence a func- 
tion of time determined by boundary conditions. Here 
the behavior of zc can be found by demanding the trans- 
formation (|18fl to be consistent with the requirement of 
a vanishing shift as discussed in more detail in ■ 

By making use of l)16|l . however, it is possible to de- 
scribe in the late time limit the profile of the radial met- 
ric component near the throat since the latter approaches 
the value rc Um = 3Af/2 there. 

Furthermore, according to (|14fl for fixed time at in- 
finity the peak in the profiles of the extrinsic curva- 
ture components arises at the throat. For Kg — —K^/2 



in the limit C 



its value there is obtained as 



.AS*,. 



2V3/9M w 0.3849/M. 



C. Origin and indicators of slice stretching 

As stated in [Tsj j , when demanding antisymmetry with 
respect to the throat maximal slices arc found where C 
being purely gauge can be chosen independently of time 
at infinity. Here the 3-mctric is given timc-indcpcndcntly 
by the initial data and the odd lapse can be written as 



OLodd = ± 



VpcF) 



(20) 
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In particular, the odd lapse vanishes at the throat and 
is positive/negative in the original/extended part of the 
Schwarzschild spacetime to yield the values plus/minus 
one at right-/left-hand spatial infinity. 

Excluding odd boundary conditions where no slice 
stretching occurs, for a discussion of the late time be- 
havior of maximal slicing instead of the slice label C it 
turns out to be convenient to introduce S as in [Tsj by 



(21) 



The late time limit t — > oo then can be said to correspond 
to the limit S — > as the maximal slices approach the 
limiting slice r = rc Hm = 3M/2 asymptotically. 

By analyzing in this limit the transformation l(18|) and 
the behavior of the 3- metric JTSJ, slice stretching effects 
can be studied. As in Q this discussion will for simplicity 
be restricted to the throat and the event horizon acting 
as markers for slice sucking and wrapping. 

In this reference it has been shown that integrating up 
to the event horizon, in the limit <5 — > the integral on 
the left-hand side of l|18|) diverges like 



V A dy 

Vpc(y) 



Cii m Cl In 



_6_ 

Jl 



0(1) (22) 



where f2 is a further fundamental constant given by 

3 



o 



'■V6M. 



(23) 



With throat and event horizon moving away from each 
other, in general also slice wrapping effects in between zq 
and z CEH show up in the form of an unbounded growth 
and/or a rapidly steepening gradient in the radial met- 
ric component. To study those, note that in numeri- 
cal implementations often a time-independent conformal 
factor \l/ 4 (z) is factored out from the 3-metric to focus 
on the dynamical features of the metric rather than the 
static singularity. In order to discuss the behavior of the 
rescaled 3-metric, it is convenient to introduce 



* 4 (z) 



which according to l|16|> are related by 



9(r, z) 



Hz) 
r 4 (r, z)' 



(26) 



Differentiating now 1)26(1 with respect to z by making 
use of the product rule and l|17|). one can extend the 
study of by analyzing in addition to the value also the 
gradient of g. In particular, at the throat one obtains 



dg 
dz 



1 dh 

r% dz 



(27) 



The proof of this statement is rather lengthy, see Q for 
details. The divergence proportional to In [S] can be un- 
derstood, however, by observing that for C — ► Cu m the 
lower limit of integration rc becomes a double counting 
root of the polynomial pc(r), the root of which appears 
in the denominator of the integrand of (|22|l . Since the 
upper limit of integration essentially plays no role in this 
expansion, one should note that any isosurface described 
by a constant value of r = const > Teh could be used 
in the following as a marker for slice stretching effects. 
Furthermore, one should note that for the constant l|23l) 
numerically a value of °f2 « 1.82M has been found in [271 
and that its analytical value has been derived in [Hi Il8j |. 

The diverging term picked up at the throat in (|22|) is 
the origin of slice stretching. Those effects are hence a 
feature of the region near the throat. Denoting the loca- 
tion of the right- and left-hand event horizon by Zq EH , 
the subscript C referring again to its time dependence, 
from the coordinate transformation (|18|l in the context 
of 1221) one can infer 



Z CEH 



± / ^H(y)dy = -C lim nki 



6 

M 



+ 0(1). (24) 



Hence slice sucking is present as in the limit S — > the 
event horizon is driven away from the throat by a term 
diverging logarithmically with 6. 



since dr/dz vanishes there. Furthermore, at the right - 
and left-hand event horizon the gradient 



dg 

dz 



'CEH 



EH 



dh 
dz 



4Ch(z± EH/ 



H{z CEH ) 



■CEH 



EH 



(28) 
is found. 

Although not obvious from the expressions (|26|l . I)27|) 
and (|28|l . one in general can expect both g and dg/dz 
to diverge in the limit of late times at the throat and/or 
the event horizon. This happens as the functions H, h 
and dh/ dz evaluated there usually grow without bounds 
while the Schwarzschild radius at the throat is approach- 
ing rc Hm = 3M/2 and at the event horizon is given by 
teh = 2M. In Scc. lIIII slicc wrapping will be worked out 
explicitly for two coordinate choices used frequently in 
numerical relativity. 

In order to describe as a function of time at infinity the 
slice stretching arising from Il22[l and showing up e.g. in 
I|24J1. it is in addition necessary to specify the relationship 
S(t) by imposing boundary conditions. Slice stretching 
for even boundary conditions will be discussed in Sub- 
sec, and favorable boundary conditions, where the 
effects described in terms of S show up late in terms of 
t, will be derived in Subsec. IIII Bl 
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III. SLICE STRETCHING EFFECTS FOR 
VANISHING SHIFT 

A. EVEN BOUNDARY CONDITIONS 



Height function 

Treating the extended and the original part of the 
Schwarzschild spacetime on equal footing by demand- 
ing symmetry with respect to the throat, even boundary 
conditions are found naturally. Here in a Carter-Penrose 
diagram the throat has to remain on the symmetry axis 
characterized by t — 0, see Fig. 4 in ||. As in by 
integrating Q) for the even "height function" one obtains 
the integral 



i(C,r) = - 



Cdy 



f(y)Vpc(y) 



(29) 



defined for r > rc- Note that the integration across the 
pole at teh is taken in the sense of the principal value 
and the corresponding slices extend smoothly through 
both the event horizon teh and the throat rc- Since 
proper time is measured at spatial infinity, from l|29(l in 
the limit r — ► oo one can infer the relationship 



Teven (C) 



between T even and C . 



Cdy 



f(y)Vpc(y) 



(30) 



e function 

The even lapse arises from (JSJ) and is given by 

dC 



± 



(34) 



dC dT even 

Studying its late time behavior as in 0] and Q , it turns 
out that a even collapses at the throat in order 0(S), 

2V2 S 



*-even \rc 



3 M 



0{8 Z ) 



0.9428— + C(<5 2 ), 



(35) 



and hence in leading order decays exponentially in time 
on the fundamental timescale VL there to avoid the sin- 
gularity. By symmetry, at both left- and right-hand 
event horizon in the limit of late times the finite value 
C Hm /r 2 EH = 3^3/16 w 0.3248 is found, 



"■even \teh 



16 



1 



7V6 In 



144^/2 
+330- 36V6- 12> 

S 2 



9^6 + 22 
3\/2 + 4 

M 2 



0{5 3 ) 



0.3248+ 1.2265 



M 2 



Q{5% 



(36) 



and by construction the even lapse approaches unity at 
both infinities. For a particular choice of spatial coordi- 
nates, namely isotropic grid coordinates, the time evolu- 
tion of the even lapse is shown in Fig. ^ 



Late time analysis 

As shown in [Tsl ] by expanding T even in terms of 5 in the 
late time limit C — > Cu m , i.e. for S — > 0, time at infinity 
is diverging like 



Teven{5) = -fHn 



S 

M 



A + 0(8) 



(31) 



where has been defined already in (|23|l and A is a 
constant with the analytic value given by 



A 

M 



V6\n [l8(3V^-4) 



2 In 



3\/3 - 5 
9\/6- 22 



(32) 



Here the divergence of T even proportional to In [5] arises 
at the throat for the same reason as pointed out for the 
expansion (|22fl . 

As later on expansions in 5 will be studied, one should 
observe that solving in (|31|l for 5, in leading order with 



6 

Jl 



exp 



cxp 



C(exp 



o 



) (33) 



an exponential decay of 5 with r e . 
timescale Q is found. 



on the fundamental 



Slice stretching in the limit of late times 

The discussion of slice stretching is particularly sim- 
ple for even boundary conditions as due to an isometry 
condition the location of the throat in Eulcrian coor- 
dinates is given time-independently by its initial value 
z\/c = z eh = const. From the expansion (|24|l in the con- 
text of O one can then infer 



C EH 



± 



y r mj)dy = C l 



%ra ' even 



O(i) 



(37) 



Hence in the limit of late times slice sucking is present 
since the event horizon is driven away from zeh by a 
term diverging proportional to r even . 

From l|26(l one can then see that very little evo- 
lution is present at the throat since g only grows 
to (4/3) 4 « 3.1605 times its initial value there as 
the Schwarzschild radius declines from teh — 2M to 
r c !im = 3M/2. Furthermore, by a symmetry argument 
the gradient of g at the throat vanishes. 

In order to study slice wrapping at the event horizon, 
it is then best to look at specific spatial coordinates such 
as logarithmic or isotropic grid coordinates which both 
are used frequently in numerical simulations. 
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First example: Logarithmic grid coordinates 



Second example: Isotropic grid coordinates 



Logarithmic grid coordinates rj [28j arise when imple- 
menting initially the Schwarzschild geometry in terms of 
logarithmic coordinates corresponding to the 3-metric 

^ds 2 = * 4 (7 ? )(d?7 2 + dfl 2 ). (38) 

= with the conformal factor given by 



Here at r e 



#(77) = \/2cosh 



'M 



the relationship 



0,77) = * 2 (r?) 



between r\ and the Schwarzschild radius 
Independently of time, the throat is 
?7vc = Veh = and the isometry 



(39) 



(40) 

r is found. 

located at 



'/ 



(41) 



is mapping the right-hand part of the spacetime to the 
left-hand part and vice versa. 

With G = '3> 12 /r 4 one may readily verify 



H(rj) = #^(77) and h{rf) = * 8 (r/) 



(42) 



and observe that the rescaled radial component of 
the metric grows from unity to the finite value 
(4/3) 4 w 3.1605 at the origin being the location of the 
throat, where in the limit of late times g behaves like 



Furthermore, since g = ^ H /r is even, obvi- 



ously its derivative vanishes there. 

Discussing slice sucking at the event horizon, from l|37|) 
and (|42|) it turns out that in leading order the event hori- 
zon moves outward like 



Vceh - ± 3 ln 



24C, 



HmTeven 



M 3 



O In 



-1/3 



(43) 



Inserting this result in (|26[) while using l|42|l . slice wrap- 
ping is taking place there as g in leading order grows 
according to 



* 8 ('/, 



CEH) 



O 



EH 



(r 4/3 ) 

\ even J 



(44) 



Furthermore, making use of l|28|) and again (|42|l . it 
turns out that a rapidly steepening gradient at the event 
horizon is present as 



dg 
dr] 



^CEH 



4Q m * 14 (r7) 



O 



EH 



(r 7 J! n ) (45) 



is found. 

These analytical statements should be compared with 
numerical results as e.g. in 0,0]. Note that in these 
simulations it is rather slice wrapping than slice sucking 
which causes the runs to crash quite early. This can now 
be understood by the argument that the event horizon 
is moving outward only moderately whereas a rapidly 
steepening gradient in g is found there together with a 
peak growing slightly further inside. 



Isotropic grid coordinates x [23 have been constructed 
in such that the 4-metric coincides at all times with 
output from a numerical evolution of black hole puncture 
data. Here initially the 3-mctric 



< 3 >ds 2 



^ 4 (x)(dx 2 +x 2 dn 2 ) 



(46) 



is implemented making use of isotropic coordinates. So 
with the conformal factor 



M 

*z =1 + — 

2x 



it turns out that x at r e 
Schwarzschild radius by 



(47) 

is related to the 



r{r e 



0,x) = x^ 2 {x). 



(48) 



Since for even boundary conditions also during the evo- 
lution isotropic and logarithmic grid coordinates are re- 
lated by 



M 



(49) 



one can see that the region 77 < is compactified to 
< x < M/2 and ?/ > is mapped to x > M/2. Here 
due to the isometry 



Ax 



(50) 



the throat is fixed for all times at xyc = x eh = M/2 and 
the puncture at x = is simply a compactified image of 
spatial infinity. 



From G = x "J 12 /r 4 then follows 

H(x) = xH 12 {x) and h(x) = .T 4 * 8 (a;) 



(51) 



and as for logarithmic grid coordinates one can observe 
that g = x 4 ^> s /r 4 at the throat grows from unity to 
(4/3) 4 w 3.1605 while due to the isometry a vanishing 
gradient is present there. 

When analyzing slice stretching at the event horizon, 
it turns out that for the outward-movement of the right- 
hand event horizon in leading order by using (|37|) and 
(I5T|) (or alternatively (g3J> an d (03) for its location 



"CEH 



(3(3 



vm 1 even 



1/3 



o 



-1/3 



(52) 



is obtained. The left-hand event horizon, however, due 
to the isometry (|50|l approaches the puncture like 



M 2 



CEH 



4.c 



CEH 



O (r- l jf) . (53) 



For both right- and left-hand event horizon the rescaled 
radial metric component g coincides and in leading order 
diverges according to 



9\x ± 



c ±4 

CEH 



f 8 (j 



ceh) 



o 



EH 



(r 4/3 ) 

\ ' even I 



(54) 
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Analyzing now the gradient of g at the event horizon, 
the particular problem of isotropic grid coordinates in 
the context of even boundary conditions comes to light. 
Whereas according to (|28|) at the right-hand event hori- 
zon the gradient 



dg_ 

dx 



+ 6 
CEH 



= o< 



r + 

'CEH 



(55) 



EH 



is found, at the left-hand event horizon the derivative 



dg_ 

dx 



dx CEH 

dx j^j dx 



-CEH 



M 1 
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is diverging even more rapidly. 

As can be seen in Fig. ^ a numerically very cumber- 
some "double peak" in the profile of g is developing which 
due to slice wrapping in the compactificd left-hand part 
of the Schwarzschild spacctime prevents long-lasting sim- 
ulations. In Subsec. 1111 HI an analysis of other bound- 
ary conditions will show, however, that initially select- 
ing isotropic coordinates can nevertheless be a "good" 
coordinate choice, since the numerically unfavorable be- 
havior described so far can be blamed mainly on the use 
of even boundary conditions. Applying for the puncture 
evolution more adapted "zgp" boundary conditions, i.e. 
demanding symmetry at the puncture and obtaining a 
vanishing gradient of the lapse there, results in signifi- 
cantly better slice stretching behavior as can be seen in 
Fig.H 

With the Schwarzschild radius at the throat approach- 
ing the value rc lim = 3M/2 in the limit of late times, by 
making use of g = :r 4, J ,8 /r 4 it is (independently of the 
boundary conditions) possible to describe the limiting 



profile of g near the throat by x 



The latter 



has been plotted in both Fig. ^ and [21 as dotted line. 

In addition, extending the study of 0], the profiles of 
the radial and the angular component of the extrinsic cur- 
vature shall be discussed here (again for arbitrary bound- 
ary conditions). According to 1)14(1 . for fixed time at infin- 
ity with r oc - for {r — > oo, x — > 0} both K r r and Kg near 
the puncture are of order 0(x 3 ), whereas with r = x for 
{r — > oo, x — > oo} for large values of x they decay in or- 
der 0(x~ 3 ). Furthermore, the peak in the profiles of the 
extrinsic curvature components is found at the throat. 
In the late time limit for Kg = —K^/2 its value there is 
obtained as C Hm /r% Hm = 2V3/9M w 0.3849/M, the lat- 
ter being in excellent agreement with numerical results 
as shown in both Fig. ^ and [2 



B. FAVORABLE BOUNDARY CONDITIONS 

Lapse constructed as a superposition of odd and even lapse 

With the trace of the extrinsic curvature vanishing, 
K = 0, the lapse arises for maximal slicing from the el- 
liptic equation 



Aa = V'V.a = Ra 



(57) 



where R is the 3-dimensional Ricci scalar. For fixed time 
at infinity, this condition is a second order linear ordi- 
nary differential equation. Hence, demanding the lapse 
to be one at spatial infinity in order to measure proper 
time there, when supplementing an additional boundary 
condition the lapse is completely determined. 

With odd and even lapse two linearly independent 
lapse functions satisfying (|57|l have been found as pointed 
out in 2 ■ By the superposition principle it is then pos- 
sible to construct a new lapse, normalized again to unity 
at right-hand infinity, by a linear combination 

a(r, r) = $(t) • cw„(r, r) + (1 - $(r)) • a odd (T, r) (58) 

with a time-dependent "multiplicator function" $(r). An 
important example for such a superposition is the punc- 
ture lapse constructed in [(|. 



Height function 

As shown in Subsec. [II Al the maximal slicing condition 
fixes the partial derivative of t with respect to r as in Q 
only, whereas boundary conditions have to be specified 
to obtain t(r, r) by integration. The latter can always 
be written as the sum of the even height function and 
a "time translation function" tc( T ) depending on time 
only, 



i(r,r) =t even (T,r) +tc(r), 



(59) 



where time at infinity is measured again in the limit 
r — > oo. As for r = one starts with the time-symmetric 
t = hypcrsurfacc, the function tc vanishes initially and 
is determined during the evolution by boundary condi- 
tions. Furthermore, since the even height function van- 
ishes at the throat during the evolution, tc also repre- 
sents the value of t at the throat. Hence the time trans- 
lation function describes where the throat is found in a 
Carter-Penrose diagram, see for "zgp" boundary condi- 
tions Fig. 4 in 0. 



Slice stretching integrals 

Imposing boundary conditions, the multiplicator func- 
tion $(t) in 1)58(1 . the time translation function tc( T ) m 
(159(1 and the location of the throat in terms of Eulcrian 
coordinates are determined. Deriving zc as a function 



8 



of t explicitly, however, is rather involved as one has to 
examine each boundary condition separately when an- 
alyzing the coordinate transformation (|18|) while mak- 
ing sure that the shift vanishes 0. For a study of slice 
stretching at e.g. the event horizon, though, the location 
of the throat has to be determined first, since, as shown 
in Subsec. Ill CI the diverging term in the integral i|24|> 
being proportional to In [6] and causing slice stretching 
is picked up at the throat. 

Integrating metric quantities from the left- to the right- 
hand event horizon yields an alternative approach for a 
discussion of slice stretching while avoiding in an elegant 
way inconveniences involved in determining zq as a func- 
tion of time. Although such integrals can not provide 
items of information like the location of throat or event 
horizon and value or gradient of the radial metric com- 
ponent there, they are nevertheless excellent indicators 
for the overall slice stretching. 

Here two such integrals shall be introduced, namely 



z + 

Z CEH 



S H (S) 
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and 



Z CEH 



S G (S) 



V G(t, y) dy 



Z CEH 
TEH 



y 2 dy 
y/Pciv) 



'C Hrl 



hi 



S 

M 
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(61) 



Whereas the first integral arises directly from the coordi- 
nate transformation 118|) and essentially has been stud- 
ied in I122[I and H24|) already, the second integral is in a 
more straightforward way related to slice sucking (since 
the left- and right-hand event horizon appear as limits 
of integration) and slice wrapping (since the integrand 
is the root of the radial metric component). So whereas 
for some analytical purposes making use of Sh might be 
preferable, for numerical studies of slice stretching Sc 
should be of particular interest. Here, in order to locate 
both left- and right-hand event horizon on the maximal 
slices to calculate either ijrJU)) or ijfiTl) numerically, one can 
compute the Schwarzschild radius r from the prefactor of 
the angular part of the metric as e.g. in 0,1201 and iden- 
tify an event horizon as isosurface with r = teh = 2M. 
Alternatively, since for the Schwarzschild spacctimc event 



and apparent horizon coincide, for this task both appar- 
ent and event horizon finders can be used in principle. 

Furthermore, one should note that the shift function 
does only enter in the calculation of these slice stretch- 
ing integrals by the fact that for zero shift the function 
H by construction depends on the spatial coordinate z 
only, whereas H (like G) is a function of both r and z 
otherwise. For the discussion of Sec. HVl it is essential to 
observe here already that the overall slice stretching as 
defined by 16U|) and l|61|l does not depend on the choice 
of the shift function. 



Boundary conditions for late observation of slice stretching 

The idea now is to obtain relationships between 6 and 
r which make the overall slice stretching, in terms of S 
arising from (|22[1 and expressed by integrals such as l|6U|) 
or (|61|) . occur late in terms of r in numerical simulations. 
By specifying S(t), however, boundary conditions for the 
lapse arise since the multiplicator function appearing in 
the linear combination (|58|) can be written as 



dS . 



dr, 



dr, 



dS 



dS 



dr e 



dr 

V dS dS J d5 dS 

This expression one can readily verify making use of |(HJ), 
(|58|l and 159fl . Assuming that a given numerical code for 
a chosen resolution can only handle a certain amount of 
overall slice stretching, the hope then is that longer last- 
ing evolutions covering a greater portion of the spacetime 
can be obtained by imposing more favorable boundary 
conditions instead of even ones. 

As one can see from (|62|) , demanding that slice stretch- 
ing effects show up later than for even boundary con- 
ditions implies with —dr/dS > —dT even /dd > that the 
multiplicator function has to be less than one, < $ < 1. 
Note that slice stretching can occur arbitrarily late with 
$ approaching zero and the lapse arising in (|58J) being 
then essentially given by the odd lapse. Furthermore, if 
$ for odd boundary conditions vanishes throughout the 
evolution, no slice stretching at all is found. 

In particular, when looking at r = T even +tc, by mak- 
ing use of l|31(l one can readily verify that if the time 
translation function tc is of order 0(1), as for the even 
lapse an exponential decay of 5 with r is found. The over- 
all slice stretching then can be expected to be similar to 
the one arising for even boundary conditions. 

Adding, however, a term which diverges logarithmi- 
cally with 8, i.e. tc — —Clin [S/M], the exponential decay 
of 5 takes place on a time scale given by the sum of S7 
and Q. The corresponding multiplicator function for this 
choice of tc is given by <& ~ S1/(S7 + ft). In principle with 
S7 — > oo the new time scale can be made arbitrarily large 
to make slice stretching effects occur arbitrarily late. In 
this limit, however, one finds that with $ — > the odd 
lapse is approached. 

It is also possible to obtain in leading order not an 
exponential but a power-law decay of 6 with r, allowing 



(62) 



9 



for very moderate slice stretching behavior. Assuming 
that t c ~ 5~ k , k > 0, then r ~ 5~ k and $ ~ flS k /k are 
found. But note that for late times with 6 — ► again the 
odd lapse is approached. 

Whereas it is possible to find boundary conditions such 
that the overall slice stretching occurs arbitrarily late, 
the corresponding lapse approaching the odd lapse might 
from the numerical point of view be disadvantageous as 
negative values of the lapse occur. For this reason a 
numerically favorable lapse in addition should be non- 
negative. Using H58(l and the formulas Ij20|l and (|34|l for 
odd and even lapse, one can show that demanding a > 
corresponds to <f> > 1/2. Hence the power-law decay of 5 
is ruled out, but choosing < f2 < il it is possible to find 
a lapse such that the time scale for the exponential decay 
of S is up to twice the one obtained for even boundary 
conditions. In particular, the non-negative lapse show- 
ing latest possible occurrence of slice stretching is at late 
times given by the average of odd and even lapse. 
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(67) 



At the right-hand event horizon, as for 

even boundary conditions, the finite value 

Cu m /r 2 EH = 3\/3/16 w 0.3248 is found in the limit 
of late times, 



First example: The puncture lapse 

It is at this point essential to note that the punc- 
ture lapse discussed in Q is precisely of this form and 
the puncture evolution of a Schwarzschild black hole is 
hence taking place in a numerically favorable manner. As 
shown in this reference, the multiplicator function of the 
"zgp" lapse is given by 



y(y-3M) dy 



2 j _ 



(63) 



(y-^y/p^iy) 



M 



Here the integral appearing in both the numerator and 
denominator of (|63|l is diverging proportional to 1 /S 2 as 
the throat in the limit of late times becomes a three- 
fold root of the denominator of the integrand. For this 
reason the "zgp" multiplicator function is of the form 



<I\~ 



1/2 + 0{8 ) and the puncture lapse, 



( T zgp, r ) 
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(64) 



being positive everywhere, arises at late times as average 
of odd and even lapse. Performing a late time analysis, 
it turns out that a zgp collapses to zero in order 0{5 2 ) at 
the puncture and the left-hand event horizon, 
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Finally, the lapse yields unity at infinity in order to mea- 
sure proper time there, see for details. The profile of 
the puncture lapse for the puncture evolution is shown 
in Fis 

In |6j the "zgp" height function is derived as 



zap 

(C,r) — L even (C,r)±(r e 



(69) 



where time is measured at right-hand spatial infinity by 

T zgp {C)=2Teven{C)-— . (70) 



Making use of the expansion (|31|l and C = Cu m + 0(S 2 ), 
in leading order the exponential decay of S with r zgp is 
with 
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taking place on twice the fundamental timcscalc Q. 

For the puncture evolution the slice stretching behav- 
ior has been discussed in from left to right at the 
left-hand event horizon, the throat and the right-hand 
event horizon. As shown there, in the limit T zgp — > oo 
the left-hand event horizon is found at a finite value in 
between its initial location xeh — M/2 and the punc- 
ture, whereas both throat and right-hand event horizon 
are moving outward like 



xc 



2 ^limTzgp 



1/3 



zgp 
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and 
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respectively. 



As 



'CEH 



freezes at late times, also the rescaled radial 



metric component approaches a finite value there. At 
xc and Xq EH , however, with g = x 4 ^> s /r 4 the metric 
diverges as in leading order 



and 
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(75) 



are found. 

Extending the study of 0, it turns out that also 
the gradient of g at the left-hand event horizon freezes, 
whereas making use of (|27|l and (|28|l the derivatives 



dg 

dx 



■re 



Ax^_ 



0(t z 



and 



dg_ 

dx 



4Clj m X CEH _ , 2 . 

„9 W Vzgp) 

'EH 



(76) 



(77) 



at the throat and the right-hand event horizon arc ob- 
tained. 

Comparing these "zgp" late time statements with the 
corresponding ones obtained for even boundary condi- 
tions in Subsec. Till Al i.e. comparing Fig. [21 with Fig. ^ 
one can see that in leading order identical slice sucking 
and wrapping is present at the right-hand event hori- 
zon and to its right. In the puncture evolution, how- 
ever, almost no slice stretching occurs to the left of the 
throat since the "zgp" lapse collapses exponentially in 
time there. For this reason numerical evolutions of black 
hole puncture data imposing the "zgp" boundary condi- 
tion are able to last significantly longer than runs forcing 
even boundary conditions. 



Second example: 
A one-parameter family of boundary conditions 

As a further example, a one-parameter family of 
boundary conditions ranging from odd to even and char- 
acterized by a constant multiplicator function in the lin- 
ear combination (|58|l . <!> = const £ [0, 1], shall be studied 
numerically in the context of isotropic grid coordinates. 
Here the lapse is determined by its timc-indcpcndcnt 
value at the puncture given by 



The elliptic equation l|57|l for the lapse has been im- 
plemented in the regularized spherically symmetric code 
described in [2j|. Using a shooting method and start- 
ing at the puncture with the value l|78|) . to be inter- 
polated there as the origin is staggered in between grid 
points, the derivative of the lapse has been determined 
such that when integrating outward a Robin boundary 
condition [30l | is satisfied. All simulations shown in this 
paper have been carried out using 30, 000 grid points for 
a resolution of Ax = 0.001A/ to place the outer boundary 
at x = 30M. 

It is worth mentioning that for negative values of the 
lapse no difficulties have been encountered numerically. 
In particular, evolving for odd boundary conditions up 
to T dd = 25 M , the deviations of lapse and metric com- 
ponents from their initial profiles have been found to be 
less than 0.1 per cent. 

For the even lapse, however, note that the isometry 
condition (|50() has not been enforced actively. In addi- 
tion one should then remember that at late times the 
left-hand event horizon according to l|53|) gets arbitrarily 
close to the puncture, whereas the even lapse approaches 



i/ r EH 



3^3/16 w 0.3248 at 
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and 



a(r,i = 0) = 2$-l6[-l,l] Vr. 



(78) 



the value Cu m/ 
is one at x = as mentioned previously. Due to the 
rapidly steepening gradient close to the puncture the 
shooting method for the even lapse failed shortly after 

7~even = 25iK/ . 

In Fig. |3 on a logarithmic scale the decay of S with 
time at infinity is shown for runs with constant values 
$ = {0, 1/8, 1/4, 1/2, 1} and for the puncture evolution 
where & Z gp — ► 1/2 holds in the limit of late times. In this 
limit analytically an exponential decay on the timescale 
f2/$ is predicted, the corresponding slopes are shown in 
addition in this figure. 

Furthermore, in Fig. 0] for the same runs the slice 
stretching integrals Sh and Sq, and (|IHTl . are plot- 
ted together with the expected late time divergence being 
proportional to time at infinity. 

As one can see from these plots, the numerical results 
are in excellent agreement with analytical predictions. In 
particular, one can see that these curves do not depend on 
the shift function, as runs demanding e.g. "zgp" bound- 
ary conditions and implementingeither elliptic gamma- 
freezing |3lT ] , minimal distortion [3, llOj or minimal strain 
[ll|,|23| shift conditions yield the same results for Sh and 
Sg as simulations using zero shift. 

In addition, note that for "zgp" boundary conditions 
as compared to the run with $ = 1 /2 more slice stretch- 
ing arises initially. This happens since the puncture 
lapse starts with unit lapse everywhere whereas a "pre- 
collapsed" lapse profile is found when using the average 
of odd and even lapse from the beginning. At late times, 
however, identical slice stretching behavior is found as 
the two curves become parallel in both Fig. |31 and 0] 

In the next section analytic arguments will be given 
on how one can cure slice stretching, analyzing in Sub- 
sec. [^^]the technique of throat excision, and discussing 
for shift functions in Subsec . IIV Bl their failure and in 



11 



Subsec. IIV CI their working mechanism for logarithmic 
and isothermal grid coordinates, respectively. 



IV. AVOIDING SLICE STRETCHING 
A. Throat excision 

As studied analytically in Subsec. Hi CI and as observed 
numerically for lo gari thmic or isotropic/isothermal grid 
coordinates in [TfllllL IT^ . Hi^ or 0, OJ, IE , respectively, 
slice stretching is arising for maximal slicing in the vicin- 
ity of the throat. Following an idea attributed to W. 
Unruh in [8^ | and excising the troublesome region of the 
hypersurfaces, long-lasting evolutions have been obtained 
numerically which do not show the previously described 
slice stretching effects, see e.g. [3ll22|. 

For an analytic understanding of such excision tech- 
niques, consider hypersurfaces starting at some inner ex- 
cision boundary, given in terms of Schwarzschild and grid 
coordinates by possibly time-dependent expressions rcex 
and zcex, respectively. Here excision shall take place in 
between the throat and the right-hand event horizon, so 
rc < rcex < r E H and z c < z Ce x < Zceh n °kl Then 
one can show that integrals of the form Sh or Sq, how- 
ever with the integration starting now at the excision 
boundary, remain finite as 
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J ^y)dy = J 



V 4 dy 

\/pc{y) 



0(1) (79) 



and 



z + 

"C EH 
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(80) 



are obtained. Here one should remember that the event 
horizon as upper integration limit in (|79fl and (|8l)|l - or 
in 1601) and l|61|) - has been chosen arbitrarily and any 
other marker lying outside of the excision boundary - or 
outside of the throat - could be used for such a calcula- 
tion. Since integrals like l|79|) and l|80|l do not diverge but 
arc finite in the limit of late times, one can now conclude 
that dynamics present in the 3-metric freeze in this limit. 
This has been observed numerically e.g. in Q when im- 
plementing horizon-locking coordinates and using several 
different shift conditions. 

In the following no use of throat excision shall be made 
when discussing in Subsccs. IIVBI and IIV CI the role of 
a non-trivial shift for evolutions using logarithmic and 
isothermal grid coordinates, respectively. 



B. Shift function and logarithmic grid coordinates 

For logarithmic grid coordinates 77 the throat is fixed 
at the origin and, treating both sides of the throat on 



equal footing by imposing even boundary conditions, the 
isometry 77 < — > —77 is present during the evolution. 

Essentially by making use of the integral Sh , i|60fl i f° r 
evolutions with vanishing shift in Subsec. IlII Al the slice 
stretching effects have been worked out at the event hori- 
zon acting as a marker. For an arbitrary shift, when 
integrating the square root of the radial metric compo- 
nent over the throat up to the event horizon, by studying 
the integral 6>g, (|61|) . one can observe that even the use 
of non-trivial shifts does not change the situation fun- 
damentally: Either slice sucking has to be present as 
the event horizon is driven away from the origin, or slice 
wrapping has to occur as the radial metric component 
blows up, or a combination of both effects is present. 

Furthermore, one should note that making use of the 
time-independent conformal factor ^(ri), and evolv- 
ing (j(t, 77) = G(T 7 r/)/fy 4 (r]) instead of G(r, 77), does not 
cure the problem. This is due to the fact that the confor- 
mal factor for logarithmic grid coordinates is finite every- 
where, as opposed to the conformal factor l|47|l which for 
puncture data diverges at the puncture, see Subsec. llV CI 
The particular form of l|39|) . however, namely its expo- 
nential growth for large values of 77, is responsible for 
the fact that numerical simulations for logarithmic grid 
coordinates are supposed to crash due to slice wrapping 
(leading to large gradients in the radial metric compo- 
nent) rather than due to slice sucking (as the outward- 
movement of the right-hand event horizon becomes ne- 
glectable at late times). 

This failure of the shift function to counteract slice 
stretching when maximally slicing a Schwarzschild black 
hole in logarithmic grid coordinates has been observed 
frequently in numerical relativity, but an analytic under- 
standing of this observation was lacking so far. In par- 
ticular, the geometrically motivated elliptic minimal dis- 
tortion shift 0, |2?| which many people in the numerical 
relativity community expected to cure the slice stretch- 
ing problem, numerically had been found to fail. It is 
instructive to cite several earlier numerical works here, 
as these references point out the difficulties encountered 
numerically when using logarithmic grid coordinates: 

D. Bernstein, D. Hobill and L. Smarr reported in [Tol l 
that "while the minimal distortion shift vector does re- 
duce a measure of the distortion [...] and removes the 
sharp gradients of the radial metric component in the re- 
gion of the event horizon, the coordinates shear is trans- 
ferred to the throat [...] and produces even larger gradi- 
ents." As a consequence, their code when using the min- 
imal distortion shift terminated at about T even = 80M as 
compared to runs with zero shift lasting for more than 
200Af. Furthermore, they realized that "part of the prob- 
lem may be blamed upon the use of a logarithmic radial 
variable" . Similarly, G. Daues in [T^ | found the minimal 
distortion shift to be "troublesome in the region near the 
throat because it leads to a large amount of slice stretch- 
ing there" . Finally, D. Bernstein in a numerical study pdj 
of the minimal distortion and various other shift condi- 
tions observed that they fail for the Schwarzschild black 
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hole, since "the minimal distortion gauge is designed to 
minimize the time derivatives of the conformal metric 
components only in a volume integral averaged sense" . 



C. Shift function and isothermal grid coordinates 

In order to make statements for evolutions taking place 
in terms of the isothermal grid coordinate x, one has to 
study the transformation of the maximal slices derived 
in the radial gauge l|13|) to the line element 



ds 2 = -A 



G 



dr 2 +2Bdrdx+Gdx 2 +r 2 dQ 2 . (81) 



Here the isothermal grid coordinate x generalizes the 
isotropic grid coordinate x of previous sections by allow- 
ing for a non-vanishing shift B. One can then observe 
that 



A(t,x) = a(T,r(r,x)), 



(82) 



B(t,x) = / 3(r ) r(r,x))|^+ 7 (T ) r(r,x))^||(83) 



and furthermore 



G(t,x) = \& (x)g(T, x) = j(t, r(r, x)) 



Or 

dx 



(84) 



hold. 

In the following, the working mechanism of the shift 
function shall be studied which makes use of the diver- 
gence of the conformal factor ^>(x) = 1 + M/2x at the 
puncture. For this study the "zgp" boundary condition 
- which in Subsec. IlII Bl has been found to be the numer- 
ically favorable boundary condition with latest possible 
occurrence of slice stretching effects - shall be chosen. 
Furthermore, a model shift should act such that during 
the evolution the rescaled radial metric function freezes 
with a profile close to one everywhere, 
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the right-hand event horizon at late times is found at 
a fixed location. Note that for this analysis the shift 
has been specified by its desired action rather than by 
demanding a particular geometric condition as in || . 
With these assumptions one can then observe from 
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that both throat and left-hand event horizon are driven 
toward the puncture since the conformal factor ty(x) di- 
verges there. At late times in leading order this move- 
ment has to take place according to 
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respectively. 

Using these results, one can now discuss the late time 
behavior of the 4-metric at five markers being puncture, 
left-hand event horizon, throat, right-hand event horizon 
and spatial infinity. 

In particular, making use of the late time statements 
of Subsec. IIII Bl it is immediately possible to sketch for 
the puncture lapse A zgp the corresponding lapse profile. 
Here a "collapsing shoulder" is obtained which, however, 
does not move outward as the location of the right-hand 
event horizon has been locked. 

The rescaled radial metric component g zgp by con- 
struction behaves nicely everywhere, and the angular 
metric component is given by the squared value of the 
Schwarzschild radius r. At the puncture, x — 0, and in- 
finitely far to the right, x — > oo, the left- and right-hand 
spatial infinities are found with r 2 behaving like x~ 2 and 
x 2 , respectively. Furthermore, as both left- and right- 
hand event horizon correspond to the value teh = 2M 
and as in between at the throat the limiting value 3M /2 is 
approached, for the angular metric component the square 
of those values is obtained at the corresponding locations. 

Finally the late time behavior of the shift shall be de- 
termined at those five markers. The calculations, how- 
ever, are rather lengthy and involve formulas such as 1041 1 . 
lO and (EH)), lO together with idJ), lEH which have 
been derived previously for the lapse a zgp , the shift f3 zgp 
and the height function t zgp . 

According to l|83() and (|84|l one can write B zgp as 



B zgp = * 2 (x)^ 



gp 



G 



Or 



(91) 



and observe that it is natural to rescale the shift by the 
conformal factor, 



B zg p{r zg p, x) — ^ [x)b zgp (r zgp , x) 



(92) 
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Demanding |jS5J, the rescaled shift b zgp is zero at both 
puncture and infinity. To obtain statements at the event 
horizon, one has to insert teh = 2M in 1|91|) and note 
that the term dr/dT zgp vanishes. The value of b zgp at 
the left- and right-hand event horizon then turns out to 
coincide with the value of the lapse there. The task of 
calculating the rescaled shift at the throat is more in- 
volved. One might actually worry that b zgp diverges there 
as rc is a root of the polynomial pc( r ): the square root 
of which appears in the denominator of (|91|) . This, how- 
ever, does not happen since the term in the brackets of 
(|9"T|) also carries a factor of \/pc( r )- In order to calculate 
drc/dr zgp , the chain rule can be applied when differen- 
tiating pc(fc) = w.r.t. T zgp , obtaining 



drc 



C 



dC 



--hp 



(2r c - 3M) dr z 



(93) 



' c 



up 



Making then in addition use of (|69[) and (|7(J|) . it finally 
turns out that the rescaled shift at the throat is given by 



y{y-3M) dy , 9 / M 

(y-^)V^) 1 c 



I 



y(y-3M) dy 



(94) 



(v-S¥-) 2 y/™(v) 



1/M 



Since here the integral appearing in both the nu- 
merator and denominator diverges in the limit 
of late times, for b zqp at the throat the value 



= \/3/6 = 0.2887 is found for r zgp — -> oo. 



There are several comments one should make here re- 
garding these results. Whereas the angular metric com- 
ponent r 2 and the rescaled shift b zgp have been found 
to behave nicely, the variables r 2 /4 ,4 (i) and B zgp usu- 
ally implemented numerically can be expected to de- 
velop pathologies close to the puncture. Those quan- 
tities should blow up at the left-hand event horizon and 
the throat, which are moving toward the puncture ac- 
cording to ijSDJl and iJSSJl, in order O (T zgp ) and O (t z ~ 2 ), 
respectively. Whereas this behavior might arise for the 
model shift only, the overall result nevertheless indicates 
that a shift has to yield "large values" close to the punc- 
ture as the diverging term of the overall slice stretching 
is "hidden" there. 



CONCLUSION AND OUTLOOK 



terms of 5 in leading order in the limit 5^0 the overall 
slice stretching has been characterized by integrals such 
as Sh and Sg- For even boundary conditions and two 
particular coordinate choices, namely logarithmic and 
isotropic grid coordinates, slice sucking and wrapping has 
been worked out explicitly in the context of vanishing 
shift. 

Searching for favorable boundary conditions, it turned 
out that slice stretching effects described in terms of S 
can show up arbitrarily late in terms of r in numerical 
simulations if the corresponding lapse approaches the odd 
lapse and hence becomes negative in the left-hand part 
of the spacetime. For numerically favorable boundary 
conditions demanding in addition the lapse to be non- 
negative, the latest possible occurrence of slice stretching 
has been found to take place for a lapse being at late 
times given by the average of odd and even lapse. The 
puncture lapse is precisely of this form and the puncture 
evolution of a Schwarzschild black hole is hence taking 
place in a numerically favorable manner. The late time 
behavior of the latter has been worked out explicitly here. 
In addition, numerical simulations have been performed 
for the puncture evolution and for a one-parameter family 
ranging from odd to even, and convincing agreement with 
analytic results has been found. 

Furthermore, analytic arguments have been given on 
how slice stretching effects can be avoided by implement- 
ing e.g. excision techniques. When using a shift, however, 
it turned out to be essential that the conformal factor has 
a coordinate singularity in order to hide the diverging 
term of the overall slice stretching there. In particular, 
the failure of a shift function for logarithmic grid coor- 
dinates and its working mechanism for isothermal grid 
coordinates have been pointed out. 

It would be interesting to extend these investigations 
by detailed studies of geometrically motivated shift con- 
ditions such as gamma-freezing and minimal distortion. 
Further numerical work could also include algebraic slic- 
ings of the "1 + log" type |2l|. The latter are used 
frequently in numerical simulations and in some regards 
mimic maximal slicing. 

Finally, as slice stretching by an intuitive argument is 
often attributed to the singularity avoiding behavior of 
the slicing, for the Schwarzschild spacetime the maximal 
slices will be compared in [3j| to geodesic slices. The 
latter arise for evolutions with unit lapse and vanishing 
shift and correspond to freely falling observers which hit 
the Schwarzschild singularity within finite time. 



Slice stretching effects have been described which show 
up when maximally slicing the extended Schwarzschild 
spacetime. Excluding odd boundary conditions where 
the static Schwarzschild metric is obtained and no slice 
stretching occurs, slice sucking and wrapping have been 
shown to arise at the throat of the maximal slices. In 
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FIG. 1: The lapse function and components of both metric 
and extrinsic curvature are shown as obtained numerically 
at times r aven = {0, 5M, 10M, 15M, 20M, 25M}. The simula- 
tion uses even boundary conditions and isotropic grid coor- 
dinates corresponding to a vanishing shift function. The lo- 
cation of the throat and the left- or right-hand event horizon 
is denoted by boxes and down- or upward pointing triangles, 
respectively. In addition as dotted lines limiting curves are 
shown which hold in the limit r even — > oo near the throat. 



FIG. 2: For the puncture evolution corresponding to "zgp" 
boundary conditions, isotropic grid coordinates and zero shift, 
the geometric quantities as of Fig. are shown in time 



steps of Atz 



5M. For t z 



oo in a region near the 



throat the puncture lapse collapses to zero, the rescaled 
radial metric component can be described by x 4, & 8 /rc l . , 
the Schwarzschild radius approaches the value rc lim = 3M/2 
and the angular extrinsic curvature component has the limit 
C llm /r 3 c . = 2V3/9M « 0.3849/M. 
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FIG. 3: The numerically observed difference in between the 
Schwarzschild radius at the throat rc and its limiting value 
rc lim = 3M/2 is plotted on a logarithmic scale as a func- 
tion of time. Here the solid lines correspond to odd ($ = 0), 
"zgp" ($zgp — > 1/2 in the limit of late times) and even ($ = 1) 
boundary conditions. The three dotted curves from top to 
bottom are characterized by a constant multiplicator function 
having the value 1/8, 1/4 and 1/2. In the lower left corner the 
analytically predicted late time slopes are shown. 
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FIG. 4: The integrals Sh and Sg characterizing the overall 
slice stretching are shown for the puncture evolution and five 
members of the one-parameter family of boundary conditions 
having constant <f>. These curves do not depend on the shift 
function as shown here explicitly for runs satisfying "zgp" 
boundary conditions and using an elliptic gamma-freezing, 
minimal distortion or minimal strain shift, see text for details. 



